Interplay of topological order and spin glassiness in the toric code under random magnetic fields 
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We analyze the toric code model in the presence of quenched disorder, which is introduced via different types 
of random magnetic fields. In general, close to a quantum phase transition between a spin polarized phase and 
a topologically ordered one, we find that increasing the amount of disorder favors the topological phase. For 
some realizations of disorder, topological order can be robust against arbitrarily strong magnetic fields. In the 
case of the toric code in a random ±h field, we show that the system exhibits a quantum phase transition to 
a spin glass phase in an appropriate dual variables description. The survival of topological order in the spin 
glass phase is directly related to the percolation properties of the rigid lattice in the Edwards-Anderson bimodal 
spin glass model. According to recent numerical results for this model [Phys. Rev. B 82, 214401 (2010)], it is 
likely that the rigid lattice does not percolate and, as a result, a new intermediate quantum phase appears in the 
random-field toric code. In this intermediate quantum phase, topological order coexists with spin glassiness. 



I. INTRODUCTION 

Quantum systems can be dramatically influenced by disor- 
der. In phenomena where Anderson localization 1 takes place, 
a weak random potential can completely localize a band of 
states that is otherwise conducting. For two-dimensional sys- 
tems, in particular, accepted wisdom suggests that quantum 
phases are destroyed if disorder is strong enough 2 . 

In the present work we revisit this basic intuition in a spin 
lattice system exhibiting topological order—. Our main goal 
is to investigate how a nonlocal topologically ordered phase 
responds to local quenched disorder. To this end, we focus 
on an exactly solvable model that supports topological order, 
namely, the toric code^— . 

We first examine a particular case of variational tensor 
product approach to the ground state of the system 7 , dubbed 
wavefunction deformation^. In this case the ground state is an- 
alytically tractable and topological properties can be mapped 
out via the topological entropy 9 '- . All ground state equal- 
time correlators are equivalent to those of a classical 2 + 
dimensional Ising model and, in the uniform-field case, there 
is a continuous transition between a Z2 topological phase and 
a spin polarized phase through a critical point that falls in the 
two-dimensional (2D) Ising universality clas a 11,12 . 

In the second part we examine a Hamiltonian deformation, 
which is closer to physical intuition in that a random magnetic 
field is added to the toric code. This situation is analyzed 
via the conventional correspondence to a 2 + 1 dimensional 
classical model, whose uniform-field critical point falls in the 
3D Ising universality class. 

We find that the addition of a random magnetic field to 
a uniform one favors in general the highly entangled phase 
with topological order over the spin polarized phase (with no 
long-range entanglement). As a result, we find that if a uni- 
form magnetic field has driven the system into a spin polar- 
ized phase, switching on some disorder leads to a recovery of 
topological order. This is an instance where the addition of 
local disorder increases the genuine many -body entanglement 
in the system. 



In the case of a wavefunction deformation without a uni- 
form field component, topological order survives for any 
amount of disorder. This is not so for a Hamiltonian de- 
formation, where an infinitely strong random field induces a 
non-topological polarized phase with spins pointing along the 
local field directions. However we argue that, in this case, 
an intermediate phase emerges and we establish an intriguing 
connection between the topological properties of this phase 
and the percolation properties of the so-called rigid latticed in 
the Edwards-Anderson bimodal (EAB) spin glass modelii. In 
particular, we show that if the rigid lattice of the 2D EAB spin 
glass model does not percolate, as conjectured recently in a 
preliminary numerical investigation 1 -, then, in the intermedi- 
ate quantum phase, spin glass properties coexist with topolog- 
ical order. Due to the interplay and possible merging of topo- 
logical and spin glass properties in this intermediate phase, we 
shall refer to it as a topological spin glass. 



In what follows we briefly review the toric code model and 
its basic properties (Sec. [TTJ> and then proceed to the main parts 
of the paper, in which we analyze the influence of random 
perturbations on the topological order of the toric code. In 
Sec. |III]we examine a wavefunction deformation, which leads 
to a fine-tuned model with a tractable ground state subspace. 
In this case, we are able to calculate the block entanglement 
and hence the topological entropy of the randomly perturbed 
toric code. In Sec. |IV]we examine a Hamiltonian deformation 
and establish the connection between percolation properties 
of the rigid lattice and topological properties of the random- 
field toric code model, in which a topological spin glass phase 
appears if the rigid lattice does not form a giant cluster span- 
ning the system. In Sec. [V]we compare and contrast the two 
types of random perturbations to the toric code and develop a 
qualitative understanding of the stability of topological order 
against them. In the final section our results are summarized 
and some derivative open questions are discussed. 



2 



H. PRELIMINARIES: TORIC CODE MODEL AND 
TOPOLOGICAL ORDER 

Quantum states possessing topological order, such as those 
encountered in the quantum Hall effect, do not break any local 
symmetry 3 . An example of such a system is the toric code*, 
which is defined on a square lattice with periodic boundary 
conditions in both directions and spins-1/2 on its bonds. The 
Hamiltonian of the model, 



Htc — £ A s — Ab £ -Bp, 



is composed of a sum over the star operators 



i ■ 



(1) 



(2) 



which are applied on bonds meeting on a vertex s; and the 
plaquette operators 



lep 



(3) 



applied on bonds forming a square p. 

There is a fourfold degeneracy (on the torus) in the ground 
state corresponding to different topological sectors. If we 
denote by W\ and W2 the incontractible Wilson loops along 
the vertical and horizontal directions, respectively, the ground 
state manifold is spanned by 



|*«) = {u;X|*tc}; i,3 =0,1} 



where 



I* 



TC, 



= igi- 1/2 



gee 



3|0> 



(4) 



(5) 



Here g E G is a string of star operators, and |0) is a reference 
state in which all spins are pointing in the +Z direction. The 
star operators form an Abelian group G of order \G\ = 2 L _1 
for an L x L lattice C. The ground states are locally indis- 
tinguishable (see, e.g., Ref.— ), by which we mean that the 
reduced density matrix, p& = Tr^l^^) (vf y ; |, is the same for 
all ground states as long as the distinguished region A does 
not form a spanning cluster in the lattice and it is smaller than 
its complement A. 



III. WAVEFUNCTION DEFORMATION 

We add a fine-tuned random perturbation to the toric code, 
which is controlled by a set of parameters {/3j }, j = 1, . . . , N. 
In this case the Hamiltonian is given by 



the perturbation term is equivalent to a local magnetic field. 
The ground state can be obtained as in Refs . 11 ' 12 , where the 
clean-case equivalent was studied (fij = (3, V j E C). It 
was shown there that the model has a quantum critical point 
at /3 C ~ 0.4406868, and topological order is lost for /3 > (3 C . 

Without detailing the derivation^ (as it is not relevant 
for our purposes here), the ground state of the disordered 
model © is explicitly given by 



l*o> 



" 1/2 £ exp 

3GG 



<7|0), 



(7) 



where 



£exp 

gee 



(8) 



A. Mapping to a Classical Random-Bond Ising Model 

A generic configuration g |0) is uniquely determined by the 
set of star operators acting on the reference state |0), modulo 
the action of the product of all g G G, which is equal to the 
identity. Therefore, there is a l-to-2 mapping between G and 
the configuration space = {6} of an Ising model with de- 
grees of freedom 6 S living on the sites s of the square lattice. 
We change variables from the bonds (a spins) of the lattice to 
the vertices (stars or 6 spins) by choosing 6 S = —1 (+1) to 
mean that the corresponding star operator A s is (is not) 'par- 
ticipating' in the relevant g € G. Since each er-spin can only 
be flipped by its two neighboring #-spins, we have aj = 9 s 9 s i, 
where j labels the bond between two neighboring sites (s, s') . 
Hence, each random parameter /L is uniquely labeled by the 
pair s, s', and we can substitute j3j —> f3 ss r. 

With this change of variables, the ground state of H can be 
expressed as 



w = e4 cx p 




where 



see 



»') 




<7(0)|O>, (9) 



(10) 



The second quantity is, in fact, the partition function of a 
classical random-bond Ising model (RBIM) on a 2D square 
lattice at temperature T, with reduced nearest-neighbor cou- 
plings fi ss i = J ss > /T. As a result, all equal-time correlation 
functions that can be expressed in terms of 6 S variables are 
exactly the same in the two systems. 



(6) 



where the first sum runs over all stars and the second sum 
runs over the four spins in each star. For small values of /3j 



B. Topological Entropy and Phase Diagram 

We can now calculate the entanglement and topological en- 
tropy^, which serves as an order parameter for phase transi- 
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FIG. 1: The four bipartitions (Ai,Bi) used to compute the 
topological entropy according to Levin and Wen^, Stopo = 
lim r ,K_ ) . o[-5 , (^i) + S(A 2 ) + S(As) - S{Ai)}, where S(A) = 
—Tt(pa In pa) is the von Neumann entropy of subsystem A. 



tions into and out of topological orde r 1 V 6 ^ . Via the topolog- 
ical entropy the phase diagram of model © can be obtained 
exactly by identifying the correspondence between topologi- 
cal order and spin polarized phases with the paramagnetic and 
ferromagnetic phases of the mapped classical model, as ex- 
plained now. 

Using the bipartition scheme of Levin and Weni£, depicted 
in Fig.Q] a calculation following closely the steps of Ref. n 
leads to the topological entropy 



.5", 



topo 



= lim 

r,R— too 2L 



■'log 2 B( ff ) (11) 



gee 



where 



z?(g) 



Zd, twisted / \ 



Z2(9) 



. twisted 



(<?) 



z!(g)z!( g ) 



(12) 

Here Z^icj) represents the partition function of an Ising 
model with reduced nearest-neighbor interactions {/3 S s' }, and 
with fixed spins along the boundary of bipartitions 2 and 
3, respectively (see Fig. [TJ. Partition functions Z® 4 (g) are 
the analogs of 3 (g) for bipartitions 1 and 4, respectively, 

whereas zf ' 4 twlsted (g) differ from them in that the spins along 
one of the two boundaries are flipped. 

The sum over g in Stopo acts as a weighed average of the 
logarithmic term over all possible values of the spins at the 
boundary. Note that in Eq. ( fL2l the partitions with two bound- 
aries, and hence with non-trivial topology, are those appearing 
with two contributions (i.e., bipartitions 1 and 4). These con- 
tributions are responsible for a non-vanishing .Stopo- 

In the high-T phase of the RBIM (j3 aa i -> 0) the correla- 
tions are short ranged and the choice of boundary conditions 
affects the partition function with only exponentially small 
corrections. Therefore, in this case, 



Z 8, twisted {g)Z d {g) 



zt{g)zi{g) 



and 5 topo = 2. This corresponds to the topological phase of 
the original model ©. On the contrary, in a ferromagnetically 
ordered phase, e.g., when (3 SS > = j3 ^> 1, we have 



z?(g) » z{ 



d, twisted 



(g),Zi(g)->Z 



d, twisted 



while Z?(g)Z$(g) ~ Z$(g)Z$(g) still holdsii. This leads 
to /S^opo = 0. The behavior of 5t opo across such a transition is 



discontinuous, with a sharp jump from S, 



topo 
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FIG. 2: (Color online) Phase diagram the random-bond Ising 
mode! 14 (RBIM) with bipartite disorder as a function of the density 
of antiferromagnetic (AF) bonds p and T = 1//3. 



C. Bipartite Disorder 

Let us study a particular case of disorder in model (O, with 
probability distribution 



(l-p)<SL0..'-^) 



topo 



0. 



The ground state of the model, Eq. (0, is controlled by 
the partition function Z of a 2D RBIM with reduced cou- 
plings (3 ss i = ±1/T, where the density of antiferromagnetic 
(AF) bonds is p. The phase diagram of this model is well- 
known 19,20 and is given for completeness in Fig. [2] The solid 
line separates the paramagnetic (PM) and ferromagnetic (F) 
phases. In 2D a spin glass (SG) phase appears only in the zero 
temperature limit (i.e., for f3 ss i — > oo). It is noted that other 
salient features in the phase diagram are straightaway inher- 
ited from the physics of spin glasses. For instance, model (O 
exhibits a multicritical Nishimori point and a re-entrant behav- 
ior in the FM phase in Fig.|2](see, e.g., Ref.— , and references 
therein). 

The region of the phase diagram where the RBIM is in the 
paramagnetic phase corresponds to a topologically ordered 
phase of the quantum model, which is continuously connected 
to the unperturbed toric code. In this fine tuned model, topo- 
logical order is stabilized by disorder, and for instance at 
p = 1/2 topological order survives for all non-vanishing val- 
ues of T. The mapped system enters the spin glass phase only 
for T — > 0. However, specifically in 2D, the stiffness ex- 
ponent of the spin glass phase is negative^- and this leads to 
R = 1 —> Stopo = 2 down to T = 0, while the system 
becomes conformally invariant 2 ^. 

It is remarkable that if we prepare the system in the spin 
polarized phase, which can be done in a straightforward way 
by applying a strong uniform field j3j = f3 > /3 C , and then we 
increase the amount of disorder (that is, if we increase the den- 
sity of AF bonds p), we can induce a recovery of topological 
order in the disordered model of Eq. ©■ 
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IV. HAMILTONIAN DEFORMATION 

We now focus on the more general case of the toric code in 
a random magnetic field, 

n = H TC -Y, h ^t- (14) 

iec 

In the limit \hi\ — > oo every spin i polarizes in the local field 
direction and hence the ground state has no topological order. 
Contrary to the uniform-field case, such a strong random field 
can favor plaquettes for which B p = — 1. An upper bound to 
the stability of this fully polarized phase can be obtained by 
setting Xa = and asking how large A b needs to be to make 
it energetically favorable for some spins to point against the 
field direction so that B p — J\ iep erf — +1 everywhere on the 
lattice. Recall that negative plaquettes can only be removed in 
pairs by flipping strings of spins that join them on the dual lat- 
tice. Therefore, the energy gain for such a change, — 4As, is 
opposed by a field contribution 2 ^ where the sum is car- 
ried out along the chosen string. If the probability of hi being 
positive is the same as that of being negative, the density of 
negative plaquettes is rj = 1/2 and the typical separation be- 
tween two such defects is £ ~ 1/ y/fj, independently of system 
size. The field contribution amounts to 2£h ~ 2h/y/rj, where 
h is the average field strength (extremal field values need not 
be considered here as it is almost always possible to find a 
string that avoids them on the lattice). Once A^ > h/2r], the 
condition B p = +1 becomes enforced throughout the lattice 
and the system exits the fully polarized phase. 

A. Mapping to a Quantum Random-Bond Ising Model 

So long as As > h/2rj, the phase diagram can be studied 
by considering the low-energy effective Hilbert space given 
by the span of all a z tensor product states satisfying the con- 
dition B p = +1, as in Refi^. We can then simplify Hamil- 
tonian (fT4l by introducing operators 6* living on the sites of 
the lattice, with eigenvalues 9 z a = (— l)" s = ±1, where n s 
is the number of times that the four spins around s have been 
flipped with respect to some reference configuration, say the 
fully magnetized state (erf = +1). The value of erf for a spin 
at bond i is uniquely determined by the product of two 9 Z S on 
adjacent sites, erf = Z S Z S ,. Indeed, in the restricted Hilbert 
space, a spin can only be flipped as part of a 4-spin flip pro- 
cess around one of the two adjacent stars (otherwise it leads 
to a high energy state ~ 2 As-) If we interpret 9 Z as the z- 
component of a spin- 1/2 degree of freedom, the operator A s 
is the corresponding x-component. Therefore, in terms of the 
0-spins, the Hamiltonian H can be rewritten as 

H = -Xa J2 % - E h ^0 z J z s> , (15) 
(«,«') 

where in the second term we rename i o ss'. This is the 
transverse field RBIM on the square lattice, whose phase dia- 
gram has been studied for several instances of disorder, such 



as ±r (bipartite) bond disorder— and bond-diluted disorder—. 
The phase diagram of the first case is shown in Fig. |4] while 
that of the second is shown in Fig. [3] Each case is studied in 
turn in the following subsections. 

B. Bond Dilution Disorder and the Percolation Threshold 

Let us start for simplicity by considering the Hamilto- 
nian (fT4-b for a toric code in a random local magnetic field, 
with bond-diluted disorder according to the probability distri- 
bution 

V{hi) = (1 - p)8(hi) + vS(h % - h). (16) 

In other words, a uniform magnetic field of strength h is ap- 
plied to a fraction p of er spins chosen at random. 

Contrary to the ±r case, the dual mapping to 9 spins holds 
for all values of h, since the diluted field does not compete 
with the plaquette term B p . In the dual language of Eq. ( fT3T > 
this type of disorder leads to a bond-diluted transverse RBIM. 

The h-p phase diagram of this model is known 2 - 5 and shown 
here in Fig. [3] For bond dilution lower than a critical value 
(P > Pc = 1/2), the system exhibits a quantum phase tran- 
sition from a paramagnetic (PM) phase for small values of 
H/Xa to a ferromagnetic phase for large values of H/Xa- For 
bond dilution above the critical value, the phase transition is 
suppressed altogether and the system remains in the PM phase 
for any h. 

In the language of the original model ( TBl i, the PM phase 
maps onto a topologically ordered ground state. Therefore 
topological order survives irrespective of field strength all the 
way up to the critical dilution, i.e., for < p < p c . Beyond 
this threshold, a sufficiently large field takes the system into a 
spin polarized phase. 

This result can be understood using a percolation argument 
applied to the large h case where each spin cr; with hi ^ 
is polarized in the direction of the field. The Hilbert space of 
the system is thus effectively reduced to that of the remaining 
(1 — p)N spins. Below the percolation threshold (p c = 1/2) 
all clusters of polarized spins are finite. One can verify that the 
corresponding projection in the large h limit simply changes 
the lattice on which the toric code is defined, but does not 
alter its topological properties. From a quantum information 
perspective, this implies a high (50%) error threshold^ 2 ^. 

On the contrary, for p > p c a spanning cluster appears. 
As h is increased and the spins on the cluster polarize, the 
projected system breaks down into disconnected components 
and topological order is lost. For instance, winding loops are 
bound to intersect the percolating cluster and topological sec- 
tors can no longer be degenerate (as discussed in the case of 
bipartite disorder with ±r magnetic field, in the previous sub- 
section). 

C. Bipartite Disorder and the Topological Spin Glass 

In the ±r models, a uniform field of strength r is applied 
to all er-spins but its direction is randomly inverted for a frac- 
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FIG. 3: (Color online) Phase diagram of the random-bond Ising 
model (RBIM) with bond-diluted disorde^ (J = \a/K). 



tion p of the spins, 
distribution is 



Mathematically, the relevant probability 



V(hi) = (1 - p) 5 (hi + r) + pd (hi - r) 



(17) 



Topological order corresponds to the fixed point J = 
AaA ~ > oo (or, equivalently, to the 6> T -polarized PM phase). 
It is stable for all values of AF-bond density p, provided the 
disordered field strength obeys r<C^. 

The case of uniform field (i.e., the line p = 1) yields the 
2D transverse field Ising model, whose phase diagram is con- 
trolled by the 3D classical Ising mode& The ferromagneti- 
cally ordered phase, which corresponds to a non-topological 
phase in model ( TBI , is controlled by the J = 0, p = 1 fixed 
point, and it is stable in a finite portion of the phase diagram. 

There is a third region in the phase diagram in Fig. [4] con- 
trolled by the J = 0, p = 1/2 fixed point. It occurs when 
the field strength r is large with respect to the transverse field 
Aa (i.e., J is small) and yet a high density of AF bonds p pre- 
vents the onset of ferromagnetic order. In this case the system 
enters a quantum spin glass phased!. 

Contrary to model (O, the fate of topological order across 
this glass transition is difficult to determine via a calculation 
of the topological entropy of the system. Here we take an al- 
ternative route and propose that this is in fact a new phase of 
matter where topological order coexists with spin glass prop- 
erties. 

We consider for simplicity the case p = 1/2, where 
Eq. ([T5T l reduces to the transverse field bimodal Edwards- 
Anderson model. In the weak transverse field limit Aa /r <C 1 
it is customary to construct approximate ground state wave- 
functions by taking the tensor product states that minimize 
the classical EAB energy (Eq. (fT5T l with Aa = 0) and forming 
superpositions of any sets of such states that are connected by 
the action of & x s operators, s 6 C. These ground states were 
shown to involve an extensive number of tensor product states 
and they exhibit statistical properties similar to the PM state 26 
(i.e., the topologically ordered state in the dual a variables). 

The ensemble of all the bonds in the EAB model that are 
always satisfied or always frustrated in all ground state con- 
figurations is called the rigid lattice^. Clusters of 9 spins 
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FIG. 4: (Color online) Phase diagram of the transverse field RBIM 
as a function of p and J = Xa/t. A topological spin glass [in the 
original model J14H appears, indicated by the shaded region. 



connected by bonds in the rigid lattice are locked with respect 
to one another. Correspondingly, the original a spins associ- 
ated with the bonds of the rigid lattice are pinned to a fixed 
value in the ground state configurations. 

As we showed in Sec. lIVBl in the case of a diluted field, 
pinning the value of spins at random in the toric code model 
does not spoil its topological properties until we pin a perco- 
lating cluster of them. As a result, we see that the survival 
of topological order in the spin glass phase of the toric code 
model in a bimodal field is related to the percolation proper- 
ties of the rigid lattice in the Edwards-Anderson bimodal spin 
glass model. 

Very recent numerical results on the EAB model 15 suggest 
that the rigid lattice does not percolate in 2D, although the re- 
sults cannot be deemed conclusive. If they are confirmed 2 ^, 
our results suggest that topological order in model (TT~4-b sur- 
vives inside the quantum spin glass phase. This occurs for 
field strengths r stronger than A^ but below the threshold 
value 2As?7, and high AF bond densities p. As r is further 
increased beyond 2 Xb t], the system eventually exits the topo- 
logical spin glass phase and enters the field polarized phase 
discussed earlier. 



V. COMPARISON BETWEEN WAVEFUNCTION AND 
HAMILTONIAN DEFORMATIONS 

One can develop a qualitative understanding of the differ- 
ence between models © and ([Pil l by studying more closely 
the two terms they differ by, that is, ^ s cxp [— Ylies ^i a i] 
and — J^i h-icrf, respectively. We shall consider the case 
hi = ±h for convenience. 

The term in model ( TBI ) is a local magnetic field which fa- 
vors a unique tensor product state with an energy gap > 2h to 
any other state. When added to the toric code, which has a fi- 
nite gap, this term eventually dominates at large field strength 
and the ground state tends to the field-polarized state, where 
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topological order is lost. 

While the perturbation of model (0 favors precisely the 
same state, it does so in a subtler way, as it becomes evi- 
dent when we consider its action on a single star s composed 
of four spin-1/2 degrees of freedom, exp [— ~^2 ies hi<jf\ . In 
this case, the lowest energy state where each spin points in the 
local field direction has a gap to the first excited state propor- 
tional to e~ 2h — e~ 4 ' 1 . 

In contrast to model ( fT4] i. the energy difference now tends 
to zero as the field strength h is increased. (Note that this 
is not true for higher excited states, whose energy difference 
with respect to the lowest state, e 2h — e~ Ah and e 4/i — e~ 4h , 
diverges exponentially with h). As a result, there are expo- 
nentially many tensor product states within an energy band 
smaller than the gap of the unperturbed toric code irrespec- 
tively of the value of h. 

This argument illustrates why topological order cannot sur- 
vive for large values of h in ( TT4l i, while the perturbation in (|6]l 
is by far less disruptive. Nevertheless, the survival of topo- 
logical order for all values of h in © is a non-trivial result 
that cannot be straightforwardly deduced from these simple 
considerations. 



VI. DISCUSSION AND PERSPECTIVES 

We have shown that topologically ordered quantum systems 
respond to quenched disorder very differently from conven- 
tionally ordered systems 2 ^, in the following sense. In gen- 
eral, close to a transition out of a topologically ordered phase, 
which is inherently long-range entangled, the introduction of 
local disorder in the applied field stabilizes the topological 
phase. In some cases, such as for random field dilution of suf- 
ficiently high density, topological order is robust with respect 
to any strength of the external magnetic field. 

An intuitive framework in which the main results can be 
understood on a physical level may be provided by the quasi- 
particle excitations (anyons) in the system. A uniform field 
can be seen as a hopping term for otherwise static anyons, 
allowing them to propagate with some lattice momenta. For 
small fields the system has a spectral gap and retains topolog- 



ical order as the ground state is adiabatically connected with 
that of the toric code 11 . For large fields the gap closes and 
topological order is lost to anyons winding around the torus. 
It is then tempting to speculate that a disordered field takes 
on the role of either a pinning potential or a random hopping 
term, which localizes the anyons. Such a mechanism would 
indeed preserve topological order, as it denies the propagation 
of pair-created anyons around the torus. 

We also argued that quenched disorder in topologically or- 
dered systems might lead to a new quantum phase where topo- 
logical order and spin glassiness coexist. Fully understand- 
ing the complex nature of this intermediate quantum phase 
requires the combined expertise of both research areas. In 
particular, our results lead naturally to the following set of 
questions, which deserves further study. 

Spin glass transitions are usually detected via magnetic sus- 
ceptibility measurements. In the present context, these involve 
an odd number of 0-spins, i.e., they are intrinsically non-local 
in the original er-spin degrees of freedom. Are there any local 
observables in the original model that can detect the topologi- 
cal spin glass transition? What properties characterize the new 
phase? Do long relaxation timescales appear in the quasipar- 
ticle dynamics of a topological spin glass 32 ? Could the glassy 
nature of this intermediate quantum phase be exploited to en- 
gineer a topological quantum memory that is robust against 
thermal fluctuations^? 

In conclusion, new phenomena arise at the confluence of 
topological order and glassy physics, with potentially interest- 
ing implications for both areas and various useful applications 
in quantum information science. 
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